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ABSTRACT 

Context. The discovery of the Spite plateau in the abundances of ^Li for metal-poor stars led to the determination 
of an observationally deduced primordial lithium abundance. However, after the success of the Wilkinson Microwave 
Anisotropy Probe (WMAP) in determining the baryon density, Qsh^ , there was a discrepancy between observationally 
determined and theoretically determined abundances in the case of ^Li. One of the most important uncertain factors 
in the calculation of the stellar ^Li abundance is the effective temperature, Tas- 

Aims. We use sixteen metal-poor halo dwarfs to calculate new Teff values using the excitation energy method. With 
this temperature scale we then calculate new Li abundances for this group of stars in an attempt to resolve the '^Li 
discrepancy. 

Methods. Using high signal-to-noise (S/N ~ 100) spectra of 16 metal-poor halo dwarfs, obtained with the UCLES 
spectrograph on the AAT, measurements of equivalent widths from a set of unblended Fe I lines are made. These 
equivalent widths are then used to calculate new Tctt values with the use of the single line radiative transfer program 
WIDTH6, where we have constrained the gravity using either theoretical isochrones or the Hipparcos parallax, rather 
than the ionization balance. The lithium abundances of the stars are calculated with these temperatures. 
Results. The physical parameters are derived for the 16 programme stars, and two standards. These include T^s, log 
g, [Fe/H], microturbulence and ^Li abundances. A comparison between the temperature scale of this work and those 
adopted by others has been undertaken. We find good consistency with the temperatures derived from the Ha line 
by Asplund et al. (2006), but not with the hotter scale of Melendez & Ramirez (2004). We also present results of the 
investigation into whether any trends between ^Li and metallicity or temperature are present in these metal-poor stars. 

Key words. Galaxy: halo - Cosmology: early Universe - nuclear reactions - nucleosynthesis - abundances - stars: 
abundances 



1. Introduction 

1.1. The Lithium Problem 

The development of the Big Bang theory brought 
about predictions for the primordial abundances of 
the nucUdes, ^H, ^H e, ■*He, and ^Li (jOlive et al. 20001 
[Fields fc Sarkar 2006p . These values, determined from cal- 
culations of Big Bang Nucleosynthesis (BBN) for an as- 
sumed baryon to photon ratio ?7, are in at least partial 
agreement with those primordial abundances determined 
from observations of ^H, ^He, ''He, and ^Li relative to ^H. 
Historically this has given constraints on rj, or equivalently 
the baryon density fraction, ilBh^, which is standard BBN's 
one free parameter, where h is the Hubble constant in units 
of 100 km"'^ Mpc~^. Precisely determining U^h^ by other 
means can narrow down the range for the calculated pri- 
mordial abundances, and thus confirm whether the obser- 
vationally inferred primordial abundances are consistent. 
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The Wilkinson Microwave Anisotropy Probe (WMAP) 
achieved high precision in measuring Ob^^- Using mea- 
surements of the cosmic microwave background (CMB) ra- 
diation, the value ^ efe^ = 0.0227 ± 0.0006 was inferred 
( [Dunkley et al. 2008[ ). This value is in excellent agreement 
with that derived from BBN via the measured ^H/^H and 
''He/^H ratios, although the systematic uncertainties in the 
latter make ^He a poor discriminant. However, for other 
primordial isotopes the results are not in such good agree- 
ment ( jCoc fc Vangioni 2005[ their Fig. 4). The baryon den- 
sities inferred from '^He and ^Li abundances differ signif- 
icantly. ^Li has the greatest deviation. The deviation of 
the observationally inferred primordial ^Li abundance from 
that deduced from WMAP and from ^H via BBN has be- 
come known as the "lithium problem" . 

The ^He discrepancy can be accommodated within large 
uncertainties concerning the mechanisms in its production 
and destruction in stars ( |Vangioni-Flam et al. 2003) . It is 
not clear that the ^Li discrepancy can be explained as eas- 

iiy- 
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The primordial Li abunda nce was first deduced from 
observations of halo stars by Spite fc Spite (1982)| when 
they discovered the ^Li plateau, now known as the 
Spite plateau, with a mean abundance of 7i(^Li)/7i(H) 
= 1.12(±0.38) X 10"^°. Many new obsevationally deter- 
mined abundances have been published, the values of 
which can differ significantly, for example from the low 
value of nCuyrUB) = {1.23tail x W-^°){95% confidence 
limits) (Ryan et al. 2000), to a decidedly larger value of 
7i(^Li)/«(H) = 2.34 X IQ-io (Melendez & Ramirez 2004). 
In comparison the abundance calculated from WMAP and 
from 2h via BBN is nfLi)/n(H) = (4.15^1^:4^ x lO^i") 



dCoc et al. (2004)1) or (^Li)/n(H) (4.26^ 



0.73 
0.60 



10^10) 

dCyburt et al. 2001[ [Cyburt 2004| . It is obvious that the 
theoretical and observational values are not in accord with 
each other. The observed ^Li abundance is a factor of 2 to 
3 times lower than that given by BBN/ WMAP, far beyond 
the stated range of systematic uncertainties from observa- 
tional analysis or BBN calculations. 

Several possibilities have emerged in an attempt to ex- 
plain this ^Li discrepancy, but none of these, as yet, fully 
explains the problem. Broadly speaking these lie in the cat- 
egories: 

1. The discrepancy is due to unrecognised or underesti- 
mated systematic errors in the calculations of the ob- 
servationally inferred present-day Li abundance. 

2. The stars studied have destroyed some of the Li with 
which they were formed. 

3. Some Galactic Li was destroyed before these stars 
formed. 

4. The discrepancy is due to systematic uncertainties in 
the nuclear cross sections used in BBN calculation. 

5. Standard BBN has failed to accurately predict the pri- 
mordial abundances. 

A possible solution in the second category is diffu- 
sion. Several studies have been conducted on this sub- 



ject recently, e.g. Richard, Michaud & Richer (2005) and 



Korn et al. (2006) I'here it is suggested that the primor- 
dial Li is depleted in the star by diffusion, through gravita- 
tional settling, to layers in which the Li can not be detected. 
The models used in this process arc manipulated to min- 
imize destruction of Li due to nuclear burning. However, 
they do not completely manage to eliminate this destruc- 
tion. Although these predictions do lead to a possible so- 
lution to the Li problem, they are still subject to ques- 
tions and uncertainties. Using atomic diffusion alone does 
not recreate the plateau, but causes a drop in Li abun- 
dance in stars with T^s > 6000 K. It is also often hard 
to explain the smallness of the star to star scatter in ^Li 
using these processes. To solve these problems some form 
of turbulence at the bottom of the convective zone has to 
be invoked. The source of this turbulence is still uncer- 
tain, as is the amount needed, which can vary from star to 
star. The turbulence is also dependent on the temperature 
scale used. If a cooler Tcs scale is used then more turbu- 
lence is needed. However, observations of ^Li in the atmo- 
spheres of metal-poor turn off stars (Asplund et al. 2006, 
IGarcfa Perez et al. 2008|l rule out the use of a more turbu- 
lent scheme, as it would destroy the ^Li. It would therefore 
seem that the success of this solution to the lithium prob- 
lem rests on an accurate Tcff scale, and confirmation of 
whether ^Li is present on the surface of metal-poor stars. 
However, it should be noted that the presence of ^Li is 



still in question. It has been shown ( [Cayrel et al. 20d7| and 
IGarcia Perez et al. 2008^ that a reappraisal of previously 
derived ^Li abundance should be undertaken. The presence 
of ^Li also causes problems if ^Li undergoes nuclear burn- 
ing. Although the models of diffusion and turbulence are 
tuned so as to reduce this burning it is not clearly stated 
whether they eliminate it completely. As ^Li is destroyed 
at a lower temperature than ''Li, the presence of ^Li would 
rule out the possibility of any models in which ''Li is de- 
stroyed. 

Possible explanations in the third category are relatively 
new. One such mechanis m involving stellar processing is 
that of Piau et al. (2006) However, the amount of material 
which lihis requires to be cycled through stars creates other 
problems; the implications of such models require further 
investigation. 

It is of course possible that certain nuclear cross sections 
used in BBN calculations have been poorly determined. 
The effect of changing the yields of certain B BN reactions 
was recently considered by Coc et al. (2004) They found 
for example, that an increase of the 'Be(d,p)2'*He reac- 
tion by a factor of 100 would reduce the ''Li abundance 
by a factor of about 3 in the WMAP 77 range. This re- 
action has since been remeasured and precludes this so- 
lution ( [Angulo et al. 2005] ). There is also the possibility 
that systematic errors in the '^He(Q!, 7)''Be reaction are 
the cause of the ''Li discrepancy. This channel was con- 
sidered in detail in Cyburt et al. (2004). Although the ab- 
solute value of the cross section for this key reaction is 
known relatively poorly both experimentally and theoreti- 
cally, the agreement between the standard solar model and 
solar neutrino data provides additional constraints on vari- 
ati ons in this cross sectio n. Using the standard solar model 
of 'Bah call et al. (2001)[ and recent solar neutrino data 
( Ahme d* et al. 2004p . one can exclude systematic variations 
of the magnitude needed to resolve the BBN ''Li problem 
at the k, 95% confidence level ( [Cyburt et al. 2004D . 

The fifth category represents an interesting challenge, 
with several cosmological and particle physics possibili- 
ties arising which could affect BBN ( |Coc fc Vangioni 2005P . 
For example, the variation of the fine structure constant 
can induce a variation in the deuterium binding energy 
and could yield a decrease in the predicted abundance of 
''Li (jPmitriev et al. 2004[ fCoc et al. 2007|). The mod ifica- 
tion of the expansion rate during BBN (jSalati 2003P will 
also affect the light element abundances. Another possi- 
bility is that gravity is not described by general relativity 
but is instead attracted toward general relativity during 
the evolution of the Universe (jPamour fc Nordtvedt 1993a[ 
IDamour fc Nordtvedt 1993bp . B BN has been exten- 
sively studied in that scenario (|Damour fc Pichon 199"9l 
ICoc et al. 2006P . Finally, it may be that physics beyond the 
standard model is responsible for the post BBN processing 
of the light elements. One possibility recently discussed is 
that particle decay after BBN could lower the ''Li abun- 
dance and produce some ^Li as weU (Je damzik 2004^ . This 
has been investigated in the framework of the constrained 
minimal supersymmetric standard model if the lightest 
supersymmetric particle is assumed to be the gravitin o 
dFeng et al. 2004[ lElIis et al.2005 l [Hamaguchi e t al. 2007| ). 
Some models have been found which accomplish these goals 
([Jedamzik et al.20061 [Cyburt et al.2006[ ). However, all of 
these possibilities lack confirmation at present. Also the 
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other four categories of explanation listed above need to be 
ruled out in order for the last to gain significant favour. 

It is the first category of explanation that we address 
in this work. The largest of the uncertainties arises from 
the uncertain effective temperature scales for metal-poor 
stars, which we now examine in greater detail, before go- 
ing on to calculate new effective temperatures and lithium 
abundances for a sample of metal-poor main-sequence stars, 
which is the ultimate aim of this work. 



1.2. Effective Temperature (Tcs) Scale Problems 

In calculating the abundances of Li, A(LiJ3, from spec- 
tral measurements, the effective temperature is the most 
important atmospheric parameter, as ^(Li) has a high 
sensitivity to temperature: dA/dTcs^ 0.065 dex per 100 
K for halo main-sequence turnoff stars. The tempera- 
ture scale adopted by different authors varies consid- 
era bly. A comparison of effe ctive temperatures adopted 



by 



Ryan et al. (2001) 
differences for 



and 



Ramirez fc Melendez (2005) 
metallicities 



shows differences for very low metallicities ([I'e/HJ < 
-3) by as much as 500 K, with typical differences 
of ~ 200 K. Concerning the very metal-poor star 
G64-37 which features in the p resent paper, the dif- 
ference in te mperature between Ryan et al. (1999), at 



6240 K, and Melendez & Ramirez (2004) at 6775 K, is 



535 K. In these cases [Ramirez fc Melendez (2005) and 



Melendez fc Ramirez (2004) have the hotter temperature 
scale, and hence they infer higher primordial abundances 
as noted above. 

Several methods are routinely adopted in calculating 
the effective temperature. These include using the strong 
Balmer line wings (e.g. Asplund et al. 2006), photometric 
methods including the Infrared Flux Method (IRFM) (e.g. 
Ramirez fc Melendez 2005) and a spectroscopic method 
that utilizes the temperature dependence of the atomic- 
level populations, as given in local thermodynamic equilib- 
rium (LTE) by the Boltzmann equation. 

In the first method, the Balmer line wings are used as 
a temperature indicator as they have strong temperature 
sensitivities, whilst maintaining a low sensitivity to other 
physical parameters, i.e. log g ([Fuhrmann et al. 1992[) . Ha 
in particular has a high sensitivity to T^s, ^ 0.5 A 
in the width of the wings per 100 K, and is the least 
sensitive t o other parameters. Th i s technique is widely 
use d, e.g. Fuhrmann et al. (1994)[ 



Barklem et al. (2002 j[ 



and Asplund et al. (2006)[ It assunies that the wings are 
formed in LTE. 'i'his assumption makes it simpler to model, 
thus leading to more sturdy temperature scales. 

The preferred temperature scale in many studies is that 
derived by the IRFM. This method uses IR photometry and 
model fiuxes to constrain the temperature. The basic idea 
is to compare a theoretically calculated infrared flux with 
an observed infrared flux to derive an angular diameter, 
given by the equation: 



Jobs,IR 
/thoo.IR 



(1) 



where f^^^ jj^ and /^i^^jq are the observed and theoretical 
infrared fluxes respectively. This 9 can then be used, along 



with the integrated observed flux, to calculate an effective 
temperature, using the equation: 



(2) 



A{U) 



"(Li) 
"(H) 



+ 12.00 



where Fobs is the observed integrated flux at the Earth 
and (T is the Stefan-Boltzmann constant. This is an itera- 
tive process where the Tcff is used to refine the theoretical 
fluxes, thus refining the angular diameter and the effec- 
tive temperature. In essence the method comes down to 
simultaneously solving equations 1 and 2 to derive val- 
ues of 9 and Tcff- Due to the lower sensitivity of the 
infrared flux to temperature compared to shorter wave- 
lengths, when equation 1 and equation 2 are solved simulta- 
neously the calculated temperature for the star is better de- 
fined, and hence more accurate (see Blackwell et al. (1979) 
Fig.l). Mean random errors for recent work using this 
method are ~ 60 - 75 K, depending on the type of star 
anal ysed, from Ramirez fc Melendez (2005)] and 70 K 
from Alonso et al. (1999), Temperatures from this method 
have been used to calibrate other photometric indices (e.g. 
Magain (1987)] and Alonso et al. (1996)) 

The third method usually relies on the assumption of 
LTE. This implies that the Boltzmann equation can be 
utilized in the derivation of the effective temperatures. 
Boltzmann's equation contains an exponential term in 
Xi/T, where Xi is the excitation energy of the given energy 
level i, and T is the temperature of the gas. Each mea- 
sured equivalent width (EW) of an element, in most cases 
Fel, is used to determine an abundance for that element. 
If the temperature used in the calculation is too high, then 
an over-population of the more excited levels will be calcu- 
lated with the Boltzmann equation. This will lead to the 
inference that there are already more absorbers in higher 
levels, and hence it would then be calculated that fewer 
absorbers were needed to reproduce the measured EW. A 
dependence of abundance on excitation energy would then 
be seen due to these miscalculated populations. The cor- 
rect temperature will therefore be the one that nulls any 
trend in the plot of abundance against excitation energy 
(e.g. Peterson fc Carney 1979 ). 

However, all methods have, along with their strengths, 
significant deficiencies. [Barklem (2007)| has questioned 
whether the assumption that the wings of the Balmer lines 
form under conditions of LTE is reliable. This therefore 
casts doubt on the accuracy of the derived temperatures, 
with a possible temperature rise of 100 K compared to those 
derived in LTE. Along with this there are also uncertain- 
ties in the way hydroge n- hydrogen collisions, s o-called self- 
broadening, are treated. Barklem et al. (2000) found differ- 
ences in Tcff of up to 150 K between the their theory and 
the Ali-Griem tre atment of self-broadenin g. This difference 
was also found by Bonifacio et al (2007) using a Ali-Griem 
theory which has modified broadening coefficients, which 
allow for a good match to solar Balmer lines. The IRFM 
relies heavily on photometry and modelled fiuxes to derive 
the temperatures. It can therefore suffer from photometric 
errors. This method also has a dependency on model atmo- 
sphere. Finally, temperatures derived using the excitation 
energy technique currently depend on the accuracy of as- 
sumptions about LTE, about the structures of ID model 
atmospheres and the errors in gf values, damping values 
and equivalent width measurements. It has been shown 
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(jTheveniii et al. 1999l and |Asplund et al. 1999] ) that the as- 
sumptions of ID, LTE atmospheres are not always suitable 
for modelling radiative transfer of F e lines in metal-poor 

infer there is little sen- 



stars, although Gratton et al. (1999) 
sitivity to non-L'i'E. it is thus unc 
effects are. 



ear how strong these 



In previous work by some of the present authors (Ryan 
et al. 1999), temperatures were derived using a calibration 
of BVRI colours and medium resolution (AA ~ lA) spec- 
troscopy, tied principally to the Magain (1987) b-y calibra- 
tion of IRFM temperatures. In Ryan et al. (2000), an ad- 
justment of -1-120 K was made to the temperature scale to 



make it similar to th e IRFM scale of Alonso et al. (1996) 
Alonso et al. (1996) report similarly a mean difference o: 



112 K between their IRFM values and Magain's. The aim 
of the present work was to explore the use of the excitation 
energy technique to further constrain the effective tempera- 
ture scale. In other studies where excitation-energy temper- 
atures have been derived, only limited attention has been 
paid to quantifying the random and systematic uncertain- 
ties on those values. Given the importance of the tempera- 
ture scale to the lithium problem, we endeavor here to track 
the uncertainties more closely. 



1.3. ID-LTE vs. 3D/lD-Non-LTE 

Because the calculation of an effective temperature from the 
excitation-energy technique requires that the populations 
of the different energy levels be known, it is necessary to 
investigate non-LTE effects in metal-poor stars and their 
effect on the temperature scale calculation. 

As discussed by Asplund et al. (1999) the level popula- 
tions for the Fe lines are affected by departures from LTE, 
as radiation fields higher than the local Planck function 
lead to over-ionisation of the atom compared to the LTE 
case. That is, applying Saha's equation with the local value 
of the temperature at some depth in the atmosphere will 
underestimate the degree of ionisation. Moreover, different 
energy levels will be affected differently by non-LTE, so it 
is not possible to apply a single ionisation correction to 
all levels of a given ionisation state. A second factor that 
could affect excitation equilibr ium is that of th e use of 3D 
atmospheres. It can be seen in Asplund (2005) (his Figure 
8) that there can be a significant x-dependent effect on 
the abundances produced. However, in the same paper it is 
stated that non-LTE corrections appear to have the oppo- 
site effect to the 3D corrections and these corrections also 
become larger when 3D atmospheres are used. This could 
lead to 3D non-LTE becoming very similar to ID non-LTE. 
It is not known, however, that this is always the case. Non- 
LTE effects, 3D effects, and their combination relative to 
ID-LTE will therefore have to be studied to assess the exact 
impact on the temperatures derived by this technique. 



However, it is well known that non-LTE calculations for 
Fe I in metal-poor stars are quite uncertain, principally be- 
cause of uncertainty in the ro le of collisions with hy drogen. 
Thevenin & Idiart (1999) and |Gratton et al. (1999")] reached 
opposite conclusions about the importance, or otherwise, of 
non-LTE in such stars. 3D calculations are also relatively 
rare to date. Because of these uncertainties, we proceed in 
two stages. For the present paper we shall concentrate on a 
ID LTE analysis, and this will be the focus for the remain- 



der of this paper. In a second paper we will examine the 
complexities introduced by attempting a non-LTE analysis. 



2. Observations 

2.1. Sample Selection 

The sample of stars was taken from one constructed by 
Ryan et al. (1999) which contains 24 metal-poor stars 
spread around the sky. It is a carefully selected sample 
at the Population II main-sequence turnoff, with as few 
star-to-star variations as possible. This included limiting 
the effective temperature to a range of 6100 K ^ Tctf ^ 
6300 K and the metaUicity to -3.5 <, [Fe/H] <, -2.5. These 
restrictions allow the sample to have roughly the same evo- 
lutionary state, in this case around the MS turn-off, which 
also has the effect of limiting log g. The low metaUicity 
restrictions mean that we are sampling stars that contain 
material that has undergone a limited amount of processing 
since the Big Bang. The sample was previously studied to 
determine not only Li abundances, but also those of mag- 
nesium (|Arnone et al. 2005)) . based on red spectra. 

The intention in the present work was to obtain blue 
high-resolution spectra for a subset of these stars. This 
was to enable us to use the excitation-energy technique to 
determine the T^s, and determine what, if any, correction 
is required to the photometric Tes scale used previously. 
This also allows for a comparison b etween the excitation- 
ene rgy Teff and those scales used by Asplund et al. (2005) 
and|Melendez fc Ramirez (2004)" 



2.2. Data Acquisition and Reduction 

Spectra of 16 of the program stars plus two standards, 
out of the 22 program stars and two stand ards, HD140283 
and HD74000, in the 
tained using the UCL 



Ryan et al. (1999) study, were ob- 
l)S instrument on the 3.9m Anglo- 
Australian Telescope (AAT) during two separate observing 
runs. The first run, July 31-August 3 2006, was undertaken 
by S.G. Ryan and J.D. Tanner. The second run, March 9- 
11 2007, was undertaken by A. Hosford and A.E. Garcia 
Perez. The setup was the same for both observing runs: 
the 79 lines/mm grating was used, centred at 4174 A with 
a slit width of 1.0", which achieves a resolving power of 
R = A/AA = 43000. The wavelength range is from 3700 A 
- 4900 A. In Table [1] we present the basic information in- 
cluding the photometric (and partially spectroscopic) tem- 
perature, Tphot (Ryan et al. 1999) for each of the stars, 
along with the S /N per 0.025 A pixel near 4300 A from the 
present work. 

The data reduction was done using standard routines 
within the IRAF package. 



3. Data Analysis 

Analysis of the data proceeded in, broadly, two stages. 
Firstly, the equivalent widths of spectral lines of Fe were 
measured. For the second stage we use WIDTH6 (Kurucz 
& Furenlid 1978) to calculate the abundances for each line, 
from which the physical parameters can be deduced. The 
remainder of this section discusses these stages. 
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Star name 


RA (1950) 


DEC (1950) 


V 


B-V 


S/N @ 4300 


Tphot (K) 


Initial [Fe/HJ 


BD- 13" 3442 


Uh 


44m 


18s 


-13"49'54" 


10.37 


0.33 


80 


6210 


-2.73 


BD+20°2030 


08h 


13m 


13s 


19°51'24" 


11.20 


0.40 


95 


6200 


-2.64 


BD+24°1676 


07h 


27m 


39s 


24°11'42" 


10.80 


0.40 


100 


6170 


-2.38 


BD+26°2621 


14h 


52m 


00s 


25°46'12" 


11.01 


0.50 


95 


6150 


-2.54 


BD+26°3578 


19h 


30m 


29s 


26°17'06" 


9.37 


0.43 


110 


6150 


-2.24 


BD+3°740 


04h 


58m 


38s 


04° 02' 24" 


9.82 


0.38 


100 


6240 


-2.70 


BD+9°2190 


09h 


26m 


35s 


08°51'24" 


11.14 


0.36 


95 


6250 


-2.83 


CD-24° 17504 


23h 


04m 


39s 


"24°08'42" 


12.18 


0.32 


100 


6070 


-3.24 


CD-33 1173 


03h 


17m 


34s 


— OO Ui zi 


10.91 


0.39 


100 


6250 


—2.91 




21h 


30m 


48s 


~oO oy Iz 


lU.Oo 


n Q7 
U.o / 


1 hp; 
lUo 


DUDU 


— Z.OO 


OU-71 1234 


16h 


02m 


18s 


—71 14 00 


10.44 


0.46 


95 


6190 


—2.60 


G64-12 


13h 


37m 


30s 


00° 12' 54" 


11.49 


0.41 


110 


6220 


-3.24 


G64-37 


13h 


59m 


53s 


-05°24'18" 


11.14 


0.36 


100 


6240 


-3.15 


LP635-14 


20h 


24m 


13s 


-00°47'00" 


11.33 


0.47 


115 


6270 


-2.66 


LP815-43 


20h 


35m 


21s 


-20°36'30" 


10.91 


0.39 


95 


6340 


-3.00 


HD 84937 


09h 


46m 


17s 


13°59'18" 


8.28 


0.42 


110 


6160 


-2.12 


HD 140283 


15h 40m 22s 


-10°46'18" 


7.24 


0.46 


170 


5750 


-2.54 


HD 74000 


03h 38m 31s 


-16°09'36" 


9.671 


0.43 


150 


6040 


-2.02 



3.1. Equivalent Width Measurement and Data Quality 
Checks 



Equivalent widths (EW) of apparently unblended Fe I and 
Fe II lines in the spectrum of each star were measured by 
fitting Gaussian profiles using the IRAF 'splot' program. 
Graphs of the Gaussian FWHM in velocity space (i^whm) 
were plotted against the measured equivalent width of the 
line, see Figure ([1]). Due to the nature of echelle spectra, in 
which AA/A is constant throughout the spectrum (where 
A A is the width of a pixel in wavelength units), it is possible 
to use these plots as checks of the quality of the measure- 
ments. The weak lines, EW ^ 90 mA, form a plateau due 
to the constancy of the FWHM of the Doppler core of lines 
on the linear part of the curve of growth. The Doppler 
core is of course broadened by the instrumental profile, 
which likewise has a constant width in velocity units. For 
intermediate strength and stronger lines, the plot curves 
upward due to the increase in their line wing strength, 
leading to larger FWHM. Lines whose equivalent widths 
fell significantly outside of these trends were re-measured 
or discarded as blends or noise-dominated measurements, 
and were rejected from further calculations. The equivalent 
widths measured in each star are listed in Table Al, A2 
and A3 in the appendix to this paper. (The appendix is 
available online only.) 

A comparison between equivalent w idths of the Fe lines 
in the star CD-24° 17504 measured by [Norris et al. (2001) 
and remeasured by A. Hosford showed good precision. 
The standard deviation of the difference between the two 
sets of equivalent widths was 1.3 mA with a mean of 
—0.58 mA for lines measured from the same spectra. Lines 
which were measur ed in two spectra, one obtained from 
Norris et al. (2001)| at S/N « 100, and one obtained by A. 
Hosford at S/JN « 50, were also compared. This comparison 
gave a standard deviation of 5.8 mA and a mean of —0.39 
mA. To increase the accuracy of the EW measurements for 
the star CD-24° 17504 we used a combination of EW mea- 
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Fig. 1. Plot of wpwHM versus equivalent width for the re- 
tained lines in the star HD140283 



the inverse of the uncertainty in the EW (ctew)- This (Tew 
is calculated using the equation: 



CTEW 



AAp,^VM 
S/N 



(3) 



surements from both the Norris et al. (2001) spectra and 
the A. Hosford spectrum, weighting each measurement by 



where l^Xpix is the width of a pixel in wavelength units, 
M is the number of pixels comprising the full width of the 
spectral line and S/N is the signal-to-noise of the spectra. 



3.2. Abundance Calculations 

To calculate the abundances from the measured lines we 
use the single line radiative transfer program WIDTH6. 
This program calculates the EW of an unblended absorp- 
tion line for a ID model atmosphere. It varies the starting 
abundance, iteratively calculating an EW until that equals 

the measured EW for the line. 

The KURUCZ06 ( http : / /kurucz . harvard . edu /grids . html| 
downloaded in Nov 2006) model atmospheres have 
been adopted for this work as it has been stated 
(jCastelli et al. 19"97ll that these no-overshooting models 
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are preferable in deriving accurate temperatures to earlier 
Kurucz overshooting models (jKurucz 1993 1). This is due 
to the effect the different methods of treating convection 
have on the temperature structure of the atmosphere, and 
therefore on the spectral features. The no-overshooting 
models lead to a better fit to observable spectral features 
at different temperatures. These models are interpolated in 
Tcff , log g and [Fe/H] over the ranges 3500 < Tcff < 50000, 
0.0 < \og g < 5.0 and -0.5 > [Fe/H] > -4.0. A com- 
parison between temperatures derived using the different 
overshooting models can be seen in Sect. 3.3. 

The log gf values of the Fel lines were com- 
piled 



from several sources: 



Blackwellet al. (1979) 



Ijlackwell et al 
tively the 



Oxford Pel 



Blackwel 



Blackwell et al. 
Blackwell et al. 



et al. 



consortium 



Bard et al. (1991) 



Bard & Kock (1994) 
tEe 



and 



1979 



1980 
coUec 



Puhret al. (1988 



O'Brian et al. 1991 



some cases 
adopted, see 
values are 



mean 
Table Al 
compiled 



T'hevenin et al 
of several 
for details 
from: 



log gf 
The Fell 



values 



m 
is 

log gf 



Hannafo rd et al. (1992 

Heise fc Koc k (1990)1 IKroll fc Ko ck (1987)[ [Moity (1983 
and Schnabel, Kock fc Holweger (1999), Values given by 
Moity (1983)1 are corrected by —0.06 dex, for those lines 
with upper-level energy < 48000 cm~^, and —0.24 dex 
for lines with upper-level energy great er than the above 



value; this follows the suggestion by Fuhr et al. (1988) 
Da mping values h ave been calculated using the tabulations 



of jAnstee fc (TMara (1995), Barklem fc O'Mara (1997)] 
and [Barklem et ar^igggJnRy;MrXIi)98)| noted that the 
damping treatment has a significant effect on the excita- 
tion energy temperatures because of the x-dependence of 
damping treatments. 

3.3. Adopted Approach to Constraining Parameters 

Several other physical parameters of a star, in addition 
to temperature, could in principle be constrained using 
WIDTH6, given appropriate data. These are log g, metallic- 
ity and microturbulence. Metallicity is determined from the 
calculated Fe abundance. Log g could be constrained by en- 
suring that the abundances derived for Fe I and Fe II lines 
are equal, and microturbulence is found by making sure 
there is no trend between the measured equivalent widths 
of the lines and the abundance. 

It would, in principle, seem possible to simultane- 
ously constrain both log g and Toff- However, as previ- 
ously mentioned, there are concerns about the non- validity 
of the assumption of LTE for Fel lines in metal-poor 
turnoff stars. This would therefore make it dangerous to 
try to constrain Tcff at the same time as log g. This 
concern became evident when Toff and log g where si- 
multaneously constrained via a WIDTH6 analysis for the 
star HD140283, using Fel and Fell hues. For this star 
results converged at Teff/log^/pe/Hj/^ = 5573/3.1/ - 
2.67/1.5 using the KURUCZ06 modelfl However, we know 
that the Hipparcos parallax gives a reliable value for 
log g for HD140283: log g ^ 3.73 ± 0.11 for a mass 
of 0.8360 M(7) as taken from the Yonsei-Yale isochrones 



(http://csaweb.yonsei.ac.kr/'^kim/yyiso.html). The uncer- 
tainty is dominated by the parallax error. Clearly the log 



^ We also tried this with the KURUCZ93 models and found 
5439/2.8/-2.79/1.5. 



g value we have derived is too low. This star has also 
appeared in many other studies leading to temperatures 
around 5750 K (Table 2), approximately 200 K hotter than 
we derived. It can then be seen that for this star, attempt- 
ing to constrain Tes and log g simultaneously from Fei 
and Fell lines drives the gravity and possibly the tem- 
perature values down to unusually low levels. Non- LTE ef- 
fects may be responsible. It is, however, still viable to use 
WIDTH6 to constrain the Toff- This is because the main 
non-LTE correction involved in Fe is expected to be ove- 
rionisation, which may significantly affect the abundance 
of Fel and therefore have a substantial effect on log g 
when calculated through ionisation balance. The correc- 
tion on X is very roughly the same for all levels within 
a particular star. This can be seen from Figur es 2 and 3 
from [Collet et al. 20 051 while [Asplund (2005) | also states 
this view. jShchukina et al. (2005)[ do hnd that tne effect can 
be different for difl'erent excitation potentials, however, this 
seems far more pronounced when using 3D atmospheres, 
compared to using ID atmospheres. Therefore overionisa- 
tion will have less affect on the trend between [Fe/H] and 
X, and it is still acceptable to use this to constrain Tcff, 
subject to the caveats in Sect. 1.3. 



Table 2. List of results for log g and Tcff of the star 
HD140283 derived using different techniques. 



Study 



Temperature 



Tcff /logs 



This work 



As plund et al. (2006) 



Rya n et al. (1996; 



Alonso et al. (1996) 



Simultaneous fit 5573/3.1 ^ 
Balmer line wing 5753/3.7^ 
Photometry 5750/3.4 ~ 

IRFM 5691/4.0^ 



" The uncertainties are 75 K and 0.15 dex respectively 
The uncertainties are 30 K and 0.04 dex respectively 
The uncertainties are 69 K and 0.5 dex respectively 



It is therefore preferable to constrain log g indepen- 
dently of deducing T^s, rather than using Fe lines to con- 
strain both. There are several options for doing this. In 
rare cases the Hipparcos parallax could be used. However 
a large proportion of the program stars do not have these 
data available, and for the stars where it is available the 
parallax errors are large. Alternatively, WIDTH6 could be 
used with a fixed temperature, such as an assumed photo- 
metric temperature, Tphot from previous work, to constrain 
log g through ionisation balance. This method still leads to 
concerns with the non-validity of LTE and the sensitivity 
of the method to the assumed temperature. As a third op- 
tion, theoretical isochrones for old, metal-poor stars near 
the main sequence turnoff can be used. This method, how- 
ever, leads to a range of possible values for log g, mainly 
due to the uncertainty in age and evolutionary state, i.e. 
whether a given star is pre- or post-turnoff. 

Because of the various concerns raised above, we 
adopted the following two stage procedure to analyse the 
stars. The first step is to use the Hipparcos parallax of the 
star if available. Table [3| contains gravities calculated from 
the Hipparcos parallax, log ^Hippj for each star, where pos- 
sible, with masses taken from the Yonsei-Yale isochrones 
with an assumed age of 12 Gyr and a metallicity suitable 
for the star. We obtain two Hipparcos gravities due to the 
different stellar masses the star could have depending on 
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whether it is on the main sequence (MS) or subgiant branch 
(SGB). In practice, the difference between these two values 
is small compared to the uncertainty arising in the parallax, 
and so we quote at most one value of log (?Hipp each star. 
Many of the stars' parallaxes have extremely large errors, 
to the point where they would have been unsatisfactory to 
use, and are excluded. In the case of the HD stars, we have 
taken the log ^fjipp the final gravity as they have small 
errors. However, for the other five stars with log (?Hipp 
Table[31 we have taken log (/jjipp to be only an interim value. 
A theoretical gravity, log g-^^^, is then used to determine a 
temperature as described below. We have used Tphot to cal- 
culate the Hipparcos gravity, log jjjipp i but note that it has 
only a low sensitivity to Tphot of +0.04 dex per 100 K. 

We also run WIDTH6 with a fixed Tphot (see Table 1) to 
find an interim (ionization-balance) log g, which we call log 
5Tphof The results of this analysis are found in Table [31 
The sensitivity of log ^Tphot to Tphot is +0.15 dex per 
100 K. The difference between the gravities, log ffxphot " 
log ^Hipp; ranges from -0.18 to -0.52, with a mean of -0.37. 
This suggests that this method, using the Tphot and the Fe 
ionisation balance, underestimates the gravity by ^0.4 dex. 
This could be due to Tphot being too cool by ~270 K, or 
may be evidence of non-LTE effects (overionisation) affect- 
ing Fel. We make no further use of the log ^Tphot- 

Because of the lack of reliable Hipparcos parallaxes, the 
unreliability of log ^xphotJ and concerns about trying to 
constrain T^s and log g simultaneously with WIDTH6, we 
then switched to using the log g^^^ implied by the Yonsei- 
Yale theoretical isochrones. These isochrones can be seen 
in Figure [21 Three different metallicities are plotted; the 
other available metallicities are [Fe/H] = -2.3 and -2.76. 
We use Tphot to obtain log g^^^, and then null the trends of 
abundance with excitation energy to obtain the excitation 
energy temperature, T^. In most cases this leads to two 
scales for temperature: a MS value and a SGB value. Only 
if the Hipparcos gravity is sufficiently accurate can we be 
sure of the evolutionary state 

An initial comparison between results calculated using 
the KURUCZ93 models, with overshoot, and those for the 
new KURUCZ06 models, with no-overshooting, was carried 
out for the star HD140283. This was done using the calcu- 
lated Hipparcos gravity for the star, log g = 3.73. It was 
found that the KURUCZ06 models gave hotter by 41 K 
than the KURUCZ93 models. We decided that for the re- 
mainder of the analysis the KURUCZ06 models would be 
used; see Sect. 3.2 for a discussion on the reasons for this. 

4. Results 

The physical parameters for the sixteen programme stars 
and the two standards are given in Table 4. 

The effective temperatures for all but the HD stars in 
Table 4 were found using isochronal gravities, log g^^^, based 
in turn on photometric temperatures. In order to under- 
stand how reliable our T^ determinations are, it is essen- 
tial that we quantify the various sources of possible error 
that contribute. As stated above, we defer discussions of 
non-LTE effects to Paper II. One source of error is the 
adopted gravity, log g^pp, for the HD stars and log g^^^ 
for the remainder. Choice of isochrones older/younger by 
1 Gyr would change log g-^^^ by ^^0.03 dex for MS and 
^0.06 dex for SGB stars, corresponding to changes in 
by 12 K and 24 K respectively. The isochronal gravity is 



3.0 




5.0 r 1 

6800 6600 6400 6200 6000 5800 5600 5400 

T.„- 

3.0 I A 




5.0 r 1 

6800 6600 6400 6200 6000 5800 5600 5400 

Fig. 2. Theoretical isochrones used in the determination of 
log g. Different ages are represented in the top panel, with 
the dashed line 13 Gyr, the sohd line 12 Gyr, and the dot 
dashed line 11 Gyr. The bottom panel gives examples for 
different [Fe/H]. The dashed line is -3.3, the solid line -2.5 
and the dot dashed line -2.15 

also sensitive to the choice of Toff, for which we use Tphot- 
An error of +100 K in Tphot would produce an error of 
typically +0.06 dex in log g^^^ in the SGB gravity, and 
-0.06 dex in the MS gravity, though more at the turnoff. 
This in turn produces an error of typically +24 K in the 
SGB T^, and -24 K in the MS T^. Hence, and importantly, 
the inferred T^ is only weakly dependent on an imperfect 
value of Tphot- Concerning microturbulence, we see that an 
error of 0.1 kms~^ leads to, on average, an error of «60 K 
in T^. Typical errors in microturbulence for this work range 
from ~0.05 - 0.15 km s^^ 

The major contributions to the uncertainty in comes 
from the procedure of nulling the dependence of [Fe/H] on 
X- This involves making a linear least squares fit to the n 
pairs of x, [Fe/H] values for a star, using n Fe lines, as may 
be seen in Figure[3l As there is a spread of [Fe/H] values for 
a given star, the slope can only be determined to a certain 
statistical accuracy, which turns out to be of order a « 0.01 
- 0.015 dex per eV, which corresponds to a temperature un- 
certainty of order ^--^60 - 80 K. The range in [Fe/H] values 
encapsulates the impact of random (not ^-correlated) line- 
to-line errors in equivalent widths, g/, and damping values. 
We determine this value for each star, finding a particularly 
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Table 3. Parameters deduced using the photometric temperatures and the gravities calculated using the Hipparcos 
parallax. Errors in log ^xphot from the ionisation balance using WIDTH6, whilst those in log gjjipp ^'"^ from the 
errors in the Hipparcos parallax. The error in log ffxphot comes about from the statistical error in making sure the the 
abundance from different ionisation states is the same, ionisation balance. 



Star name 


No. of Fel 


Tphot (K) 


log STphot 


[Fo/H] 


C(Km/s) 


log ^Hipp 




lines measured 








BD- 13" 3442 


86 


6210 


3.77 ± 0.04 


-2.72 


1.30 




BD+20°2030 


98 


6200 


3.985 ± 0.04 


-2.62 


1.20 




BD+24°1676 


96 


6170 


3.68 ± 0.04 


-2.60 


1.40 




BD+26°2621 


107 


6150 


3.88 ± 0.10 


-2.69 


1.20 


4.28 ± 0.31 


BD+26°3578 


126 


6150 


3.72 ± 0.08 


-2.39 


1.35 


3.94 ± 0.30 


BD+3°740 


94 


6240 


3.68 ± 0.12 


-2.77 


1.50 


4.18 ± 0.21 


BD+9°2190 


146 


6250 


3.68 ± 0.10 


-2.84 


1.50 




CD-24° 17504 


69 


6070 


3.57 ± 0.30 


-3.35 


1.30 




CD-33°1173 


126 


6250 


3.87 ± 0.10 


-3.04 


1.50 


4.33 ± 0.27 


CD-35° 14849 


69 


6060 


3.92 ± 0.07 


-2.44 


1.20 


4.44 ± 0.25 


CD-71°1234 


127 


6190 


3.90 ± 0.07 


-2.54 


1.50 




G64-12 


59 


6220 


4.05 ± 0.08 


-3.40 


1.40 




G64-37 


71 


6240 


4.29 ± 0.04 


-3.24 


1.30 




LP635-14 


163 


6270 


3.90 ± 0.10 


-2.49 


1.50 




LP815-43 


157 


6340 


3.78 ± 0.12 


-2.74 


1.40 




HD 84937 


164 


6160 


3.66 ± 0.10 


-2.32 


1.30 


3.98 ± 0.12 


HD 140283 


124 


5750 


3.40 ± 0.04 


-2.53 


1.50 


3.73 ± 0.12 


HD 74000 


84 


6040 


3.77 ± 0.03 


-2.19 


1.20 


4.03 ± 0.18 


Sensitivity to +100 K in Tphot 




+0.15 








±0.04 



Table 4. The final physical parameters deduced in this work for the main sequence (MS) scale and the sub giant (SGB) 
scale. The three HD stars had gravities determined by their Hipparcos parallaxes and as such their evolutionary state 
can be deduced. 



Star name 


log fisc 

(MS) 




Tx (K) 
(MS) 


[Fe/HJ 
(MS) 


e (Km/s) 
(MS) 




log ftsc 

(SGB) 




Tx (K) 
(SGB) 


[Fe/HJ 
(SGB) 


e (Km/s) 
(SGB) 


BD-13" 3442 


4.46 


± 





14 


6321 ± 87 


-2.66 


1.20 


3 


79 


± 





14 


6186 ± 87 


-2.74 


1.30 


BD+20° 2030 


4.46 


± 





14 


6208 ±104 


-2.63 


1.00 


3 


79 


± 





14 


6099 ± 104 


-2.69 


1.20 


BD+24° 1676 


4.47 


± 





14 


6296 ± 55 


-2.53 


1.20 


3 


78 


± 





14 


6220 ± 55 


-2.57 


1.40 


BD+26° 2621 


4.49 


± 





14 


6233 ± 107 


-2.67 


1.20 


















BD+26° 3578 
















3 


76 


± 





14 


6148 ± 81 


-2.42 


1.35 


BD+3° 740 


4.44 


± 





14 


6344 ±113 


-2.72 


1.50 


3 


81 


± 





14 


6188 ± 113 


-2.82 


1.50 


BD+9° 2190 


4.44 


± 





14 


6486 ± 129 


-2.68 


1.50 


3 


81 


± 





14 


6352 ± 129 


-2.75 


1.50 


CD-24° 17504 


4.55 


± 





14 


6102 ± 232 


-3.33 


1.00 


3 


61 


± 





14 


6110 ± 232 


-3.31 


1.30 


CD-33° 1173 


4.44 


± 





14 


6386 ± 55 


-2.94 


1.50 


















CD-35° 14849 


4.52 


± 





14 


6168 ± 38 


-2.36 


1.00 


















CD-71° 1234 


4.46 


± 





14 


6194 ± 53 


-2.52 


1.20 


3 


78 


± 





14 


6172 ± 53 


-2.54 


1.50 


G64-12 


4.49 


± 





14 


6333 ± 90 


-3.32 


1.30 


3 


67 


± 





14 


6304 ± 90 


-3.31 


1.40 


G64-37 


4.48 


± 





14 


6175 ± 106 


-3.29 


1.20 


3 


69 


± 





14 


6181 ± 106 


-3.30 


1.40 


LP635-14 


4.43 


± 





14 


6319 ± 114 


-2.48 


1.50 


3 


83 


± 





14 


6135 ± 114 


-2.61 


1.50 


LP815-43 


4.40 


± 





14 


6529 ± 107 


-2.61 


1.40 


3 


87 


± 





14 


6400 ± 107 


-2.68 


1.40 


HD84937 
















3 


98 


± 





12 


6168 ± 102 


-2.34 


1.30 


HD 140283 
HD74000 


4.03 ± 0.18 


6070 ± 127 


-2.20 


1.20 


3.73 ± 0.12 


5769 ± 39 


-2.54 


1.50 



large value for CD-24°17504 on account of the small num- 
ber of Fe I lines in this very metal-poor star. The final la 
uncertainty in , given in Table [31 includes the quadra- 
ture sum of the (dominant) statistical error arising from 
nulling the [Fe/H],x trend, and contributions for Aage = 
1 Gyr, = 0.1 kms"\ A[Fe/H] = 0.05, and an allowance 
of A Tphot = 100 K. 

We comment briefly on the larger scale of the uncer- 
tainties quoted here, compared to the relatively small val- 
ues (30 - 40 K) quoted for Tphot by |Ryan et al. (1999)[ The 
aim of the previous work was to derive temperatures that 
minimised the impact of possible systematic errors. The 
purpose was to investigate the spread of Li abundances by 
ensuring random errors were kept to a minimum, while at 



the same time acknowledging that large zero-point errors 
might persist but be of approximatly the same size in all 
stars. In the current work, however, the emphasis has been 
to constrain the zero-point, albeit at the expense of larger 
random errors for individual stars. 

These temperatures have been used in combination with 



the Li equivalent widths measured by iRyan et al. (1999) 
and a grid of equivalent width versus abundance for dif 



ferent Toff's taken from Ryan et al. (1996)] to calculate 
new Li abundances. The grid was constructed by deter- 
mining different ^Li equivalent widths for a number of 
different Toff and [Fe/H] values. This was done by pro- 
ducing synthetic spect ra of the ^Li region with an LTE 
code originating from Cottrell & Norris (1978) The log 
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Star name 


EW (Li) 


Teff (K) (MS) 


A(U) (MS) 


Toff (K) (SGB) 


A{U) (SGB) 


BD-13" 3442 


21.0 


6321 ± 87 


2.20 


± 


0.057 


6186 ± 87 


2.10 


± 


0.057 


BD+20° 2030 


20.5 


6208 ±104 


2.11 


± 


0.068 


6099 ±104 


2.03 


± 


0.068 


BD+24° 1676 


21.1 


6296 ± 55 


2.18 


± 


0.009 


6220 ± 55 


2.13 


± 


0.009 


BD+26° 2621 


22.5 


6233 ± 107 


2.17 


± 


0.070 










BD+26° 3578 


24.6 










6148 ± 81 


2.15 


± 


0.053 


BD+3° 740 


19.5 


6344 ±113 


2.18 


± 


0.074 


6188 ± 113 


2.07 


± 


0.074 


BD+9° 2190 


14.6 


6486 ±129 


2.14 


± 


0.084 


6352 ± 129 


2.05 


± 


0.084 


CD-24° 17504 


18.1 


6102 ± 232 


1.98 


± 


0.151 


6110 ± 232 


1.98 


± 


0.151 


CD-33° 1173 


17.2 


6386 ± 55 


2.15 


± 


0.036 










CD-35° 14849 


28.8 


6168 ± 38 


2.24 


± 


0.025 










CD-71° 1234 


25.9 


6194 ± 53 


2.21 


± 


0.035 


6172 ± 53 


2.19 


± 


0.035 


G64-12 


21.2 


6333 ± 90 


2.21 


± 


0.059 


6304 ± 90 


2.19 


± 


0.059 


G64-37 


18.2 


6175 ± 106 


2.03 


± 


0.070 


6181 ± 106 


2.03 


± 


0.070 


LP635-14 


20.2 


6319 ± 114 


2.18 


± 


0.074 


6135 ± 114 


2.05 


± 


0.074 


LP815-43 


16.1 


6529 ± 107 


2.21 


± 


0.069 


6400 ± 107 


2.13 


± 


0.069 


HD84937 


24.9 










6168 ± 102 


2.17 


± 


0.066 


HD140283 
HD74000 


47.9 
22.1 


6070 ± 127 


2.05 ± 0.083 


5769 ± 39 


2.21 ± 0.025 



> 2 
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showing the trend 



Fig. 3. Plot of X versus log(^j^ 
nulling involved in constraining the Tcff 



gf values and w avelengths of the four components of ^Li 
were taken from Andersen, Gustafsson, & Lambert (1984) 
and the m odel atmospheres were interpolated from the 
Bell (1981) grid. The details of the process can be found 
in Norris et al. (1994) Table [5] presents these values. 



5. Discussion 

5.1. Comparison with other temperature scales 

With the determination of the effective temperatures using 
the excitation method we can now compare our tempera- 
ture scale to that of others, in particular to the IRFM of 



Melendez fc Ramirez (2004) the Balmer line wing method 
of Asplund et al. (2006), and the principally photometric 



temperatures of Ryan et al. (1999) 



6800 




5800 



6000 6200 6400 6600 
Trff (K) this work 



6800 



Fig. 4. Comparison of the Melendez & Ramirez (2004) 



temperatures and the temperatures from this work. Filled 
circles represent SGB temperatures and open circles repre- 
sent MS temperatures. Those stars of unknown evolution- 
ary state have two points on the plot. 



In Figure [H the compar ison with the IRFM Toff of 
Melendez & Ramirez (2004) is plotted. That scale is clearly 
hotter than ours, by an average of ^ 250 K. This is not to 
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say that either scale is right, but it does lead us to the con- 
clusion that our temperature scale will not reconcile the Li 
problem; the Melendez & Ramirez (2004) IRFM scale was 
on the borderline of doing so. Our scale will give a lower 
mean Li abundance than will the IRFM. The difference in 
temperatures could have several possible explanations. It 
could be due to problems with our use of LTE in the exci- 
tation method, causing it to calculate lower temperatures. 
An incorrect bolometric flux, fboi, calibration used in the 
IRFM could lead to that method having higher temper- 
atures. It has also been noted (Alonso et al. 1996) that 
the errors in determining the absolute IR flux calibration 
have different effects on the derived IRFM temperatures, 
depending on what photometric band is used. The effect 
of these different errors on the derived temperatures is to 
move the zero point of th e temperature scale. This cou ld be 
another reason that the Melendez & Ramirez (2004) tem- 
peratures are hotter. 



6800 



o 
o 



6600 



Z 6400 



< 6200- 



i4 



6000 



5800 




5800 



6000 6200 6400 6600 
T,ff (K) this work 



6800 



Fig. 5. Comparison of the Asplund et al. (2006) tempera- 
tures and those from this work. I'llled circles represent SGB 
temperatures and open circles represent MS temperatures. 
Those stars of unknown evolutionary state have two points 
on the plot. 

The comparison between the temperatures of 



Asplund et al. (2006) and this work is plotted in Figure 
5) Here we see that the two scales are comparable. As 
discussed in Sect. 1.2 both of these methods may suffer 
from LTE effects. A departure from LTE of only a few 
percent would increase Balmer wing temperatures by of 
order 100 K (Barklem 2007). It is speculated that the 
use of LTE in the treatment of the Fe atom may be 
inaccurate, e.g. overionisation of the element relative to 
LTE may occur (jThevenin et al. 1 999). If true this would 
also lead to the excitation temperatures changing, as they 
are based on an LTE treatment. As yet the value of the 
temperature change is not determined and the sensitivity 
of each method to non-LTE may be very different. This is 
an aspect that we intend to look into in future work. It is 
clearly of interest to know whether non-LTE corrections 



to these scales could increase both to be akin to that 
of the IRFM scale, and push the mean Li abundances 
toward that of the WMAP inferred abundances. This 
agreement between the two scales also gives an indirect 
comparison with the Balmer wing temperature scale of 
Bonifacio et al (2007) They found that their temperature 



scale was essentially the same as Asplund et al. (2006) 
which is in turn similar to ours. 
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Fig. 6. Comparison of the Ryan et al. (1999) temperatures 
and those of this work. Filled circles represent SGB tem- 
peratures and open circles MS temperatures. Those stars 
of unknown evolutionary state have two points on the plot. 



Finally, Figure [S] presents the comparison between the 



photometric temperatures of Ryan et al. (1999) and of 



this work. At lower temperatures, 7^ < 6250 K, the two 
scales are in excellent agreement. For the MS scale we have 
seven stars cooler than 6250 K. The difference, — Tphot, 
ranges from -126 K to 56 K, with a mean of -19 K. For the 
SGB scales we have eleven stars cooler than 6250 K, with 
the difference ranging from -99 K to 47 K and a mean of 
-18 K. However, at higher values, > 6250 K, we see that 
the excitation scale becomes hotter. In this case, for eight 
MS scale stars we have a range of 87 K to 262 K and a 
mean of 143 K, and for three SGB scale stars a range of 
57 K to 117 K with a mean of 73K. This suggests that the 
Tphot values are reliable for T"^ < 6250 K, but may be ~70 
- 140 K too cool for ^ 6250 K. Photometric temper- 
ature errors can be induced through inaccurate deredden- 
ing, inaccurate calibration to standards, or a flawed colour- 
temperature calibration. The last of these is perhaps the 
most readily encountered; earlier in this work we cited dif- 
ferences of ^100 K between th e IRFM scales of Magain 
(1987) and Alonso et al. (1996) Of course, the difference 
could also be due to inaccuracies in the model atmospheres 
used in this work to derive TV. 
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5.2. Lithium Abundances 

Concerning the lithium problem, with our temperature 
scales we achieve a mean Li abundance of ^(Li)= 2.16 
dex for the MS scale and ^(Li)= 2.10 dex for the SGB 
scale, with a scatter of 0.074 dex and 0.068 dex respec- 
tively. The mean Li abundance determined from the five 
stars whose evolutionary state is defined is ^(Li)= 2.18 
dex with a scatter of 0.038 dex. It is pleasing to note 
the similari ty between our mea n ''Li abundances and that 
derived by [Spite et al. (1996)[ ^(Li)= 2.08 (±0.08) dex, 



who also use a temper ature scale based on iron excitation. 
Bonifacio et al (2007) also found a similar result, ^(Li)= 
2.10 (±0.09) dex, with their Balmer line wing temperature 
scale. 

These values are obviously still too low to reconcile 
with the WMAP value of ^(Li)«2.62 dex. Even if we later 
find that non-LTE corrections give temperature rises of 
roughly 150 K, we would get mean lithium abundances of 
^(Li)«2.26 dex and ^(Li)«2.20 dex for the MS and SGB 
scales respectively. This value is still too low. However, it 
is not yet known what is the true effect of non-LTE and/or 
3D corrections on our temperature scale. 

Although the primary purpose of this analysis was to 
investigate possibly large systematic errors in temperature, 
we can also investigate whether there is evidence for depen- 
dencies of ^(Li) on Tcs and [Fe/H] within the sample. With 
the ^(Li) values from Table 5 and the physical parameters 
calculated for each star we are able to perform a multiple 
regression fit, linear in T^s and [Fe/H]. In these fits we in- 
clude the stars with known evolutionary state, three on the 
MS and two on the SGB, combined with the eleven remain- 
ing stars of unknown evolutionary state for which we assess 
two cases: all MS and all SGB. This produces the equation 

^(Li) = 1.36(±0.84)+0.00019(±0.00013)Toff 
+0.138(±0. 048) [Fe/H] 

for the MS case and 

^(Li)=0.79(±1.22)+0.00025(±0.00019)Tcff 
+0.095(±0. 055) [Fe/H] 

for the SGB case. Here the bracketed values are 1 a errors 
in the coefficients. The large errors seen in the intercept 
coefficient comes from the large extrapolation to T = 
and [Fe/H] = 0. We see no statistically significant trend in 
either T^b or in [Fe/H] for the SGB as each coefficient has 
errors of similar size to the coefficient itself. Although there 
is the possibility of a trend existing in [Fe/H] for the MS 
case (the coefficient is ~3cr) there is no certainty in this 
as it is not known whether all the stars belong in the MS 
evolutionary state. 

The lack of significant trends here, in contrast 
to the significant [Fe/H] dependence in the st udy by 
Ryan et al. (1996) (and by Asplund et al. (2006)), comes 



about because the random uncertainties on the 'iV values 



us- 
the 



are much larger than achieved by Ryan et al. (1999) 
ing photometry. We recall that the primary purpose of 
present work was to search for potentially large systematic 
errors in the temperatures, and for that purpose the larger 
random errors here are tolerable. However, the larger ran- 
dom errors greatly undermine the second use of the dataset. 

We have also performed single parameter regression fits 
in Tcs and [Fe/H]. Figure [7] presents the results of the fits 



to temperatures. The values of the temperature coefficients 
imply 

^(Li)=0.51(±1.00)+0.00026(±0.00016)rcff 

for the MS case, and 

^(Li)=0.71(±1.33)+0.00022(±0.00021)reff 

for the SGB case. There is clearly no statistically significant 
trend with temperature in either the MS or the SGB result, 
again since the temperature coefficient errors are of similar 
size to the coefficients themselves. For the fit to [Fe/H] we 
obtained the equation 

^(Li)=2.58(±0.14)+0.152(±0.049)[Fc/H] 

for the MS case, and 

^(Li)=2.35(±0.16)+0.089(±0.056)[Fe/H] 

for the SGB case. Figure [5] shows a plot of these fits. Here 
we confirm t he similarity in va lues of the metallicity coef- 

and Asplund et al. (2006), 
re SGB case do not usefully 



ficient with iRyan et al. (1999) 
but again the large errors for t 
constrain a metallicity trend. The metallicity coefficient for 
the MS case is > 3cr, but we caution again that we cannot 
be sure that all stars are on the MS. 

We have also used the fitting form as described by Ryan 
et al. (2000) to deduce the primordial Li abundance. This 
fit is of the form: 



Li/H = a +b Fe/FeQ 



(4) 



where a is the intercept and measures the primordial 
abundance directly. From this we obtain a value of ^Li/H 
= (1.18 ± 0.10) X 10-1" for the MS case and ^Li/H = 
(1.10 ± 0.1) X 10-1" for the SGB case. Comparing this 
to the value deduced from WMAP via BBN, '^Li/H = 
4.15to;45 X 10^1", we once again see we have not recon- 
ciled the lithium problem. 

We can see from this work that it has not been possible 
to find a solution to the lithium problem through the use of 
temperatures derived using t he excitation technique. The 

values are lower than the Melendez & Ramirez (2004) 
values, and more similar to the Ha and photometric tem- 
perature estimates. It would seem that one of the biggest 
outstanding problems is the assumption of LTE in the cal- 
culations. As has been stated in previous sections, we have 
begun our own investigation into the effects of non-LTE 
on abundance analysis, and therefore the effective temper- 
atures, of this group of stars. 



6. Conclusions 

The physical parameters of 16 program stars, and two stan- 
dard stars, have been calculated using WIDTH6. In partic- 
ular, we have derived excitation energy temperatures by 
nulling the dependence of abundance on excitation poten- 
tial for Fe I lines. We have compared our temperatures with 
those calcula ted in other works using diffe rent techniques, 
the IRFM of iMelendez fc Ramirez (2004)[ the Balmer line 
wing n iethod of Asplund et al. (2006)[ wnich is similar to 
that of'Bonifacio et al (2007) and the photometric method 
of |Ryan et al. (1999)] 

We have shown that the IRFM scale of 



Melendez & Ramirez (2004) is hotter than ours by an 
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Fig. 7. Li abundances as a function of temperature. The 
dotted line is the fit to the MS temperatures (open circles) 
and the dashed line is the fit to the SGB temperatures 
(filled circles). Those stars of unknown evolutionary state 
have two points on the plot. 
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Fig. 8. Li abundances as a function of metallicity. The dot- 
ted line is the fit to the MS metallicitics (open circles) and 
the dashed line is the fit to the SGB metallicities (filled cir- 
cles). Those stars of unknown evolutionary state have two 
points on the plot. 



average of ~250 K. This difference may be the effect of 
the Ffjoi calibration used in IRFM or it could be a shift of 
the zero-point as an effect of the errors induced by using 
different photometric bands in the IRFM calculations. On 
the other hand it may be the effects of assuming LTE on 
our temperature scale. 

It has been found that our methods have produced tem- 
peratures comparable with those derived using the Balmer 
line wing method. Both of these scales are too low to rec- 
oncile the lithium problem. It has been noted that the LTE 
assumptions of the Balmer line wing method may be lead- 
ing to temperatures that are on order of 100 K too low. The 
method used in this work may also suffer from problems 
with LTE, but the full effect of non-LTE on the excitation 
method used here is unknown. This is something we will 
research in future work. 

Finally we also see comparable results with the photo- 
metric temperature scale, although at higher temperatures. 



— 6250 K, there is a shift towards our scale becoming 
hotter. This suggests that the photometric temperatures 
are reliable up until ^6250 K. The most readily detected 
error in photometric calculations comes from the colour- 
temperature calibration and differences in values derived 
from different calibrations can be as much as ^100 K. 

With our derived temperatures we have calculated the 
Li abundances of the program stars. Due to large uncer- 
tainties, we do not see statistically compelling trends in 
^(Li) with either T^s or [Fe/H]. Due to the uncertainty 
in evolutionary state for the majority of the stars we have 
two mean lithium abundances, ^(Li) ~ 2.16 dex for the 
main sequence case and ^(Li) = 2.10 dex for the SGB case. 
Five of the stars do have defined evolutionary states. The 
mean Li abundance from these five stars is j4(Li)= 2.12 dex. 
These values are lower than the WMAP determined value 
of ^(Li)«2.62 dex. For the observationally deduced value 
to be comparable with the WMAP value our temperatures 
would need to increase by approximately 700 K. 

Acknowledgements. We thank J. D. Tanner for assisting with collect- 
ing the initial data at the AAT and the referee, P. Bonifacio, for his 
helpful comments which have improved the presentation of this work. 
The work of KAO was supported in part by DOE Grant No. DE- 
FG02-94ER-40823. AH, SCR, AEGP gratefully acknowledge support 
from the Royal Society under International Joint Project 2006/23 in- 
volving colleagues at Uppsala University. 



References 

Ahmed, S. N., ct al. 2004, Physical Review Letters, 92, 181301 

Alonso, A., Arribas, S., & Martmez-Roger, C., 1996, A&AS, 117, 227 

Alonso, A., Arribas, S., & Martinez-Roger, C., 1999, A&A, 139, 335 

Andersen, J., Gustafsson, B., & Lambert, D. L. 1984, A&A, 136, 65 

Angulo, C., et al. 2005, ApJ, 630, L105 

Anstee, S. D., & O'Mara, B. J. 1995, MNRAS, 276, 859 

Arnone, E., Ryan, S. G., Argast, D., Norris, J. E., & Beers, T. C. 

2005, A&A, 430, 507 
Asplund, M., Nordlund, E., Trampedach, R., & Stein, R. F. 1999, 

A&A, 346, L17 
Asplund, M., 2005, ARA&A, 43, 481 

Asplund, M., Nissen, P. E., Lambert, D. L., Primas, F. & Smith, V. 

V. 2005, lAU Symp. 228, 53 
Asplund, M., Lambert, D. L., Nissen, P. E., Primas, F., & Smith, V. 

V. 2006, ApJ, 644, 229 
Bahcall, J. N., Pinsonneauh, M. H., & Basu, S. 2001, ApJ, 555, 990 
Bard, A., Kock, A., & Kock, M. 1991, A&A, 248, 315 
Bard, A. & Kock, M. 1994, A&A, 282, 1014 
Barklem, P. S., & O'Mara, B. J. 1997, MNRAS, 290, 102 
Barklem, P. S., O'Mara, B. J., & Ross, J. E. 1998. MNRAS, 296, 1057 
Barklem, P. S., Piskunov, N., & O'Mara, B. J. 2000, A&A, 363, 1091 
Barklem, P. S., Stempels, H. C., AUende Prieto, C., Kochukhov, O. 

P., Piskunov, N., & O'Mara, B, J. 2002, A&A, 385, 951 
Barklem, P. S. 2007, astro-ph/0702222 B 
Boll, R. A. 1981 private communication 

Blackwell, D. E., Ibbetson, P. A., Petford, A. D., & Shallis, M. J. 

1979, MNRAS, 186, 633 

Blackwell, D. E., Petford, A. D., & Shallis, M. J., 1979, MNRAS, 186, 
657 

Blackwell, D. E., Shallis, M. J., & Selby, M. J. 1979, MNRAS. 188, 
847 

Blackwell, D. E., Petford, A. D., Shallis, M. J., & Simmons, G. J. 

1980, MNRAS, 191, 445 
Blackwell, D. E., Petford, A. D., Shallis, M. J. 

1982, MNRAS, 199, 43 
Blackwell, D. E., Petford, A. D. 
210, 595 

Bonifacio, P., Molaro, P., Sivarani, T., Cayrel, R. , Spite, M., Spite, 
F., Plez, B., Andersen, J., Barbuy, B., Beers, T. C., Depagne, E., 
Hill, v., Frangois, P., Nordstrom, & Primas, F. 2007, A&A, 462, 
851 

CasteUi, F., Gratton, R. G., & Kurucz, R L. 1997, A&A, 318, 841 



Simmons, G. J. 



Simmons, G. J., 1982, MNRAS, 



A. Hosford et al.: Lithium abundances of halo dwarfs based on excitation temperature. I. LTE 



13 



Cayrel, R., StefTen, M., Chand, H., Bonifacio, P., Spite, M., Spite, F., 
Petitjean, P., Ludwig, H. G., & Caffau, E., 2007, A&A, 473, 37 

Coc, A., Nunes, N. J., Olive, K. A., Uzan, J.-P., & Vangioni, E. 2007, 
Phys. Rev. D, 76, 023511 

Coc, A., Olive, K. A., Uzan, J.-P., & Vangioni, E. 2006, Phys. Rev. D, 
73, 083525 

Coc, A., Vangioni-Flam, E., Descouvemont, P., Adahchour, A. & 

Angulo, C, 2004, ApJ, 600, 544 
Coc, A., & Vangioni, E. 2005, lAU Symp, 228, 13 
Collet, R., Asplund, M., & Thevenin, F. 2005, A&A, 442, 643 
Cottrell, P. L., & Norris, J. E. 1978, ApJ, 221, 893 
Cyburt, R. H. 2004, Phys. Rev. D, 70, 023505 

Cyburt, R. H., Ellis, J., Fields, B. D., Olive, K. A., & Spanos, V. C. 

2006, Journal of Cosmology and Astro-Particle Physics, 11, 14 
Cyburt, R. H., Fields, B. D., & Olive, K. A. 2001, New Astronomy, 

6, 215 

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2004, Phys. Rev. D, 69, 
123519 

Damour, T., & Nordtvedt, K. 1993, Physical Review Letters, 70, 2217 
Damour, T., & Nordtvedt, K. 1993, Phys. Rev. D, 48, 3436 
Damour, T., & Pichon, B. 1999, Phys. Rev. D, 59, 123502 
Dmitriev, V. F., Flambaum, V. V., & Webb, J. K. 2004, Phys. Rev. D, 
69, 063506 

Dunkley, J., et al. 2008, ArXiv e-prints, 803. iarXiv:0803.0586 l 
Ellis, J., Olive, K. A., & Vangioni, E. 2005, Physics Letters B, 619, 
30 

Feng, J. L., Su, S., & Takayama, F. 2004, Phys. Rev. D, 70, 075019 
Fields, B.F., & Sarkar, S., 2006 in Yao, W.-M., & et al. 2006, Journal 

of Physics G Nuclear Physics, 33, 1 
Fuhr, J. R., Martin, G. A., & Wiese, W. L. 1988, Journal of Physical 

and Chemical Reference Data, 17, 4 
Fuhrmann, K., Axer, M., & Gehren, T. 1992, A&A, 271, 451 
Fuhrmann, K., Axer, M., & Gehren, T. 1994, A&A, 285, 585 
Garcia Perez, A. E., Inoue, S., Aoki, W., & Ryan, S. G. 2008, 

Precision Spectroscopy in Astrophysics, Proceedings of the 

ESO/Lisbon/Aveiro Conference, 9 
Gratton, R. G., Carretta, E., Eriksson, K., & Gustafsson, B. 1999, 

A&A, 350, 955 

Hamaguchi, K., Hatsuda, T., Kamimura, M., Kino, Y., & Yanagida, 

T. T. 2007, Phys. Lett. B, 650, 268 
Hannaford, P., Lowe, R. M., Grevesse, N. & Noels, A. 1992, A&A, 

259, 301 

Heise, C. & Kock, M. 1990, A&A, 230, 244 
Jedamzik, K. 2004, Phys. Rev. D, 70, 063524 

Jedamzik, K., Choi, K.-Y., Roszkowski, L., & Ruiz de Austri, R. 2006, 
Journal of Cosmology and Astro-Particle Physics, 7, 7 

Korn, A. J., Grundahl, F., Richard, O., Barklem, P. S., Mashonkina, 
L., Collet, R., Piskunov, N., & Gustafsson B. 2006, Nature, 442, 
657 

KroU, S. & Kock, M, 1987, A&AS, 67, 225 

Kurucz, R. L. 1993, lAU CoUoq. 138: Peculiar versus Normal 
Phenomena in A-type and Related Stars, Astronomical Society 
of the Pacific Conference Series, 44, 87 

Magain, P. 1987, A&A, 181, 323 

Melendez, J., & Ramirez, L 2004, ApJ, 615, L33 

Moity, J. 1983, A&AS, 52, 37 

Norris, J. E., Ryan, S. G., & Stringfellow, G. S. 1994, ApJ, 423, 386 
Norris, J. E., Ryan, S. G., & Beers, T. C. 2001, ApJ, 561, 1034 
O'Brian, T. R., WickUffe, M. E., Lawler, J. E., Whaling, W., & Brault, 

J. W. 1991, JOptSocAmB, 8, 1185 
Olive, K. A., Steigman, G., & Walker, T. P. 2000, Phys. Rep., 333, 

389 

Peterson, R., Carney, B. W. 1979, ApJ, 231, 762 

Piau, L., Beers, T. C, Balsara, D. S., Sivarani, T., Truran, J. W., & 

Ferguson, J. W. 2006, ApJ, 653, 300 
Richard, O., Michaud, G., & Richer, J. 2005, ApJ, 619, 538 
Ramirez, L, & Melendez, J. 2005, ApJ, 626, 446 

Ryan, S. G. Beers, T. C, Deliyannis, C. P., Thorburn, J. A. 1996, 
ApJ, 458, 543 

Ryan, S. G., Norris, J. E., & Beers, T. C. 1996, ApJ, 471, 254 
Ryan, S. G. 1998, A&A, 331, 1051 

Ryan, S. G., Norris, J. E., & Beers, T. C. 1999, ApJ, 523, 654 
Ryan, S. G., Beers, T. C, Olive, K. A., Fields, B. D. & Norris, J. E. 

2000, ApJ, 530, L57 
Ryan, S. G., Kajino, T., Beers, T. C, Suzuki, T. K., Romano, D., 

Matteucci, F., & Rosolankova, K. 2001, ApJ, 549, 55 
Salati, P. 2003, Phys. Lett. B,571. 121 

Schnabel, R., Kock, M. & Holweger, H. 1999, A&A, 342, 610 



Spite, F. & Spite, M. 1982, A&A, 115, 357 

Spite, M., Frangois, P., Nissen, P. E., & Spite, F. 1996, A&A, 307, 
172 

Shchukina, N. G., Trujillo Bueno, J., & Asplund, M. 2005, ApJ, 618, 
939 

Thevenin, F. 1989, A&AS, 77, 137 
Thevenin, F. & Idiart, T. P. 1999, ApJ, 521, 753 
Vidal, C. R., Cooper, J., & Smith, E. W. 1973, ApJS, 25, 37 
Vangioni-Flam, E., Olive, K A., Fields, B. D. & Casse, M. 2003 ApJ. 
585, 611 



